Half-filled Kondo lattice on the honeycomb lattice 
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The ground-state phase diagram of the Kondo lattice model on the honeycomb lattice at half- 
filling is carefully studied by using an extended mean-field theory. By treating magnetic interaction 
and Kondo screening on an equal footing, it is found that besides a trivial discontinuous first-order 
quantum phase transition between well-defined Kondo insulator and antiferromagnetic insulating 
state, there can exist a wide coexistent region with both Kondo screening and antiferromagnetic 
orders in the intermediate coupling regime. The appearance of such wide coexistence region is 
attributed to the unique linear density of state on the honeycomb lattice at half-filling in contrast to 
the case of square lattice. It is expected such a coexistence region could be realized by experiments 
of ultra-cold atoms on the honeycomb optical lattices and may be found by quantum Monte Carlo 
simulation in future since it is free of the 'minus sign'. Additionally, the stability of Kondo insulator 
to the trivial decoupled state requires a minimum strength of the Kondo coupling, which is also 
the result of the magic linear density of state. We hope the present work may be helpful for 
further studies on the interplay between conduction electrons and the densely localized spins for the 
honeycomb lattice. 



I. INTRODUCTION 

It is still a challenge to understand the emergent quan- 
tum phases and corresponding quantum criticality in 
heavy fermion systems To attack this challenging 
problem, the so-called Kondo lattice model is introduced, 
which is believed to capture the nature of interplay be- 
tween Kondo screening and the magnetic interaction, 
namely, the Ruderman-Kittel-Kasuya-Yosida exchange 
interaction, mediated by conduction electrons among lo- 
calized spins The former effect favors a nonmagnetic 
spin singlet state in strong coupling limit while the lat- 
ter tends to stabilize usual magnetic ordered states in 
weak coupling limit. There seems to exist a quantum 
phase transition or even a coexistence regime between 
these two kinds of well-defined states^ - — however, a 
more radical critical phase may have been observed in 
YbRh2(Sio. 96 Geo.o6hm 

The nature of these mentioned phenomena has been 
a long standing but controversial issue since the work 
of Doniachi^i It is noted that an extended mean-field 
theory treating the magnetic interaction and the Kondo 
screening on an equal footing, has predicted a coexis- 
tence regime of the disorder Kondo singlet and the an- 
tiferromagnetic ordered state with small staggered mag- 
netization and partially screened local moments in the 
intermediate coupling. 1 - This coexistence has been con- 
firmed by sophisticated quantum Monte Carlo (QMC) 
simulation without notorious 'minus sign' problemi 16 i 17 
Therefore, it is fair to say that such a mean-field theory 
could provide reliable physical results in the intermediate 
and strong Kondo coupling regimes when the unconven- 
tional quantum phases, e.g., quantum spin liquids, are 
irrelevant in the related problems . 22 ' 23 

Recently, much attention has been focused in strongly 



correlated physics on the honeycomb lattice since the 
low energy excitations are described by relativistic 
Dirac fermions around distinct Dirac points rather than 
usual non-relativistic Landau quasiparticle near Fermi 
surface 2^— 

So far, most of the theoretical studies in this active 
field focus on Hubbard model or its derivative with spin- 
orbit coupling, the Kane-Mele-Hubbardi 2 ^^! 2 ^^ 
How the conduction electrons interplay with the densely 
localized spins on the honeycomb lattice, which could be 
encoded in terms of the Kondo lattice model in principle, 
is largely an open problem^ 

In this work, we try to uncover the ground-state phase 
diagram of the Kondo lattice model on the honeycomb 
lattice at half-filling using the mentioned extended mean- 
field theory for an anisotropic Kondo lattice model. Our 
main result is systematically summarized in Fig. [TJ In 
weak Kondo coupling an antiferromagnetic ordered in- 
sulating state appears, and in strong coupling limit the 
Kondo insulator is ground state. In the intermediate 
regime, a wide coexistence regime is found. In some re- 
spect, our work may be considered as the first attempt 
in this direction before more sophisticated numerical ap- 
proaches can be used. 



II. THE KONDO LATTICE MODEL ON THE 
HONEYCOMB LATTICE AND MEAN-FIELD 
DECOUPLING 



The model we considered is the anisotropic Kondo 
lattice model defined on the honeycomb lattice at half- 
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FIG. 1. A possible ground-state phase diagram of the Kondo 
lattice model on the honeycomb lattice at half-filling. AFM 
refers the antiferromagnetic ordered insulating state in weak 
Kondo coupling while Kondo insulator appears in strong cou- 
pling. In the intermediate regime, a coexistence region is 
found with both Kondo screening and antiferromagnetic or- 
ders. 
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H = H t + H\\+H u 



(cj di 



2e ik ^ 2 cos V3fc y /2 and A, 



c iBi c iBi) = 2m c and -V = 

4 t Ci t ) with f{k) = e~ ik * - 
B representing two nonequivalent sublattices. It can be 
seen that m dl m c corresponds to magnetization of local 
spins and conduction electrons, respectively, while non- 
vanishing V denotes the onset of Kondo screening effect. 
Besides, since we are considering a half-filled lattice, the 
local constraint has been safely neglected at the present 
mean-field level with chemical potential setting to zero^ 
Diagonalizing the above mean-field Hamiltonian, the 
four quasiparticle bands are obtained as 



E ±± (k) = ±-j=y/F 1 {k)±F 2 (k) 



(2) 



with FiO) = J^(m 2 d + m 2 c )/4+ JlV 2 /2 + t 2 \f(k)\ 2 and 
Fa(fc) = [F^k) 2 - A(Jlm 2 m 2 J16 + J\m d m c J\V 2 /8 + 

JiV 4 /16) - J 2 m 2 c t 2 \f(k)\ 2 ]i. Therefore, the wanted 
ground-state energy at half-filling is given by 



E„ 



(E^(k) + E_ + (k)) + E . 



(3) 



h± = -y Y^( c h c n4i d it + 4i c *t4t d *i)> (!) 

i 

where H t describes conduction electrons hopping be- 
tween nearest-neighbor sites and the pseudofermion rep- 
resentation for local spins has been utilized as S" = 
h So-cr' 4a T acr' d i<y' with T a being usual Pauli matrix and 
a local constraint 4\ d i\ + 4i d U ~ 1 enforced in every 
site. Hn denotes the magnetic instability due to the po- 
larization of conduction electrons by local spins while H± 
describes the local Kondo screening effect resulting from 
spin-flip scattering process of conduction electrons by lo- 
cal moments. The Kondo screening effect has also been 
formulated in the 1 /N expansion but magnetic instability 
is difficult to treat in such a framework. 37 

Now, using the mean-field decoupling introduced by 
Zhang and Yv>^ for the longitudinal and transverse in- 
teraction term £/y, H±, respectively, it is straightforward 
to obtain a mean-field Hamiltonian 

H MF =H t +Hf F +H? F + Eo t 

H t = -t^2(f(k)cl Aa c kBcr + f*(k)cl Ba c kAa ), 

ka 
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( C kAa d kA<j + c\ Ba d kBa + h.C.), 



E = N a (2J\\m d m c + J±V 2 ), 

where we have defined several mean-field parameters as 
(d\ At d lAt - d\ Ai d iAX ) = 2m d , (<4jt<W ~ 4 B i d iBl) = 
-2m d , {c\ M c lA \ - 4ai c ^m) = -2rn c , (cJ st c iSt - 



III. THE MEAN-FIELD SOLUTION FOR 
GROUND-STATE 

Firstly, we proceed to discuss two simple but physi- 
cally interesting limits for Kondo coupling Jm and J±, 
which correspond to the antiferromagnetic ordered state 
(J|| ^> Jl) and Kondo insulating state (J|| <C J±), 
respectively.— 

For the case with Jy 3> Jj_, in general, one expects the 
antiferromagnetic ordered state to be the stable ground- 
state of Kondo lattice model on the honeycomb lattice 
due to its bipartite featured To study the possible an- 
tiferromagnetic ordered state, making use of the ground- 
state energy (Eq(3|) with assuming no Kondo screening 
existing (V = 0) and the quasiparticle spectrum (EqJ2]), 
we can easily derive ground-state energy of the antifer- 
romagnetic ordered state per site as 
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and two self-consistent equations from minimizing E y 
with respect to magnetization m d and m c , respectively. 

J\\m c (2m d - 1) = 0, 

2J \\ m - - ^rrii^ 2 + j ii m d) - J \\ m d] = o 



2A 2 



where we have used a simplified linear density of state 
(DOS) p(e) = |£|/A 2 when transforming the summation 
over momentum k into integral on energy s with A being 
high-energy cutoff.— Thus, t\f(k)\ can be replaced by |e| 
to simplify corresponding calculations. 

From these two equations, one obtains m d = 1/2 
and m c = J||[(4A 2 + J 2 /4) - J||/2]/8. Meanwhile, 
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the low-lying quasiparticle excitations in the antifer- 
romagnetic ordered state has the energy E± + (k) — 

±^t 2 \f(k)\ 2 + J|/16 and £±_(Jfc) = ±J\\m c /2 with an 

apparent gap Jy/4 around the Dirac points where f(k) — 
0. Thus, we conclude the antiferromagnetic ordered state 
we obtained is indeed an insulating state with fully po- 
larized local spins (rrid = 1/2) while conduction electrons 
only partially polarize (m c < 1/2). This feature is simi- 
lar to the previous study on square lattice, thus confirms 
the validity of our current treatment ~ 

Another interesting case appears when Jii <C Jj_. It 
is natural to expect that a Kondo insulating state arises 
in this situation for half- filling . 13 ' 15 Following the same 
methology of treating antiferromagnetic insulating state, 
we can get the ground-state energy per site for the ex- 
pected Kondo insulating state with V ^ but no mag- 
netic orders wid = m c = 



j^Kondo j t/2 [7 \ 2 

E 9 = J ^-^ A 



j2y 2)3/2 _ j3y3 



Minimizing E^ ondo with respect to Kondo hybridiza- 
tion parameter V, we obtain V = 1 — A 2 /(4Jj) which 
implies a critical coupling = 1/2 corresponding to 
vanishing V. It is noted that V oc (J± — J]_) in con- 
trast to usual mean- field result j3 = 1/2, and this can be 
attributed to the low-energy linear DOS of conduction 
electrons on the honeycomb lattice at half-filling. Similar 
critical behavior for onset of Kondo screening on the hon- 
eycomb lattice has been obtained in the study of Kondo 
breakdown mechanism as well.— As a matter of fact, the 
existence of the critical coupling J£ = 1/2 results from 
the competition between the Kondo insulating state and 
the trivial decoupled state where V = = m c — 0, its 
ground-state energy E° = — 4A/3 comes solely from free 
conduction electrons. Comparing E® and E^ ondo , one 
clearly recovers the critical coupling which justifies 
the above simple picture. 

However, since the Kondo insulating state is unsta- 
ble to the decoupled state when J± < = 1/2, one 
may wonder whether a trivial decoupled state appears be- 
tween the Kondo insulating state and the antiferromag- 
netic insulating state. It is easy to see the ground-state 
energy of the antiferromagnetic state E^ FM is always 
lower than E® for any positive coupling Jy . Therefore, 
a trivial disordered state in intermediate coupling seems 
unfavorable based on our current mean-field treatment. 

Instead, we may consider a possible quantum critical 
point (QCP) between Kondo insulating state and the an- 
tiferromagnetic state and its position can be determined 
by comparing the ground-state energies of Ef FM and 



The position of the QCP is readily obtained as J c = 
0.88A, numerically. However, in fact, one can check that 
this putative QCP is a first-order quantum phase tran- 
sition point when comparing the first-order derivative of 
jjjKondo anc j Eg FM with respect to the Kondo coupling 
J. Thus, we do not expect radical critical behaviors near 
such first-order quantum phase transition point in the 
spirit of Landau-Ginzburg-Wlison paradigm* 1 . Addition- 
ally, it is noted that such first-order quantum phase tran- 
sition has also been obtained on the square lattice^ and 
we suspect this feature may be generic for conventional 
mean-field treatment of Kondo lattice models according 
to standard Landau-Ginzburg phase transition theory^ 
However, we point out that it is a subtle issue to com- 
pare the results of a first-order quantum transition with 
numerical simulations, particularly when the first-order 
transition is a weak one^ 

Generically, in the intermediate coupling regime, a 
possible coexistence region of the Kondo insulating 
state and the antiferromagnetic ordered state cannot be 
excluded^ For searching such possibility, we have to 
resort to the full formulism of ground-state energy E g 
(Eq|3]) and three self-consistent equations derived from 
-jjrr- = 75— - = -g— - — 0. Since the possible coexistence 
region has two boundaries, which correspond to the on- 
set of Kondo screening (V ^= 0) in antiferromagnetic in- 
sulating state and antiferromagnetic order (m c ,m f / =/= 0) 
arising in Kondo insulator, we can derive two equations 
for these two distinct boundaries. 

The first boundary which means the antiferromagnetic 
order {m c ,m,i ^ 0) arising in Kondo insulator can be 
obtained by solving the following equation 



J(l 



4 J 2 ' 3A 2 

2 A 2 
-— J 3 (l - —) 3 
3A 2 1 4J 2j 



(A 2 - J 2 



A 4 
16 J 2 



h/A 2 + j 2 (i 



A 2 . 
4J 2 ' 



(5) 



jgKondo p Qr physically interesting case with J± = J\ 
J, we have 



where we have utilized the expression for the Kondo hy- 
bridization parameter V = 1 — A 2 /(4J^) and assum- 
ing the isotropic condition (Jj_ = Ju = J). Solving 
Eq|5] numerically, one finds the first boundary located 
in J c i = 0.78A which is smaller than first-order quantum 
phase transition point (J c — 0.88A), thus the expected 
coexistent scenario is stable to the previous discontinuous 
first-order quantum phase transition. 

For the second boundary, which corresponds to the 
onset of Kondo screening (V ^ 0) in antiferromagnetic 
state, one can also derive similar equation like Eq.([5]) to 
determine the location of the boundary. 

2A 2 f A 

— = J dee[P(e) + Q(e) + (P(e) - Q(e))W(e)],(6) 



6A 2 



[(4A 2 + J 2 /4)5 - (J 2 /4)^ 



-1 



J 

A 2 



(J- 



4J 2 ' 



(4) 



where we have defined three auxiliary functions P(s) = 
l/yje 2 + J 2 /16, Q(e) = 2/Jm c and W(e) = [e 2 + 
J 2 /16 + J 2 (m 2 -m c )/4]/[e 2 + J 2 (l/4-m 2 )/4] with m c = 
J||[(4A 2 + J, 2 /4)- J||/2]/8. Then, one obtains the location 
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of the second boundary J C 2 = 1.06 A by numerically solv- 
ing Eq.©. It is noted that J c2 = 1.06A > J c = 0.88A, 
which implies the system indeed tends to form a coex- 
istent state with Kondo screening and antiferromagnetic 
orders. 

Therefore, we may conclude that there may exist a 
coexistence regime (J c i = 0.78A < J < J c2 = 1.06A) in 
the intermediate coupling regime based on the discussion 
in two phase boundaries for the isotropic case ( Jj_ = Jii = 
J). When comparing to the square lattice, we find the 
range of the coexistence regime on the honeycomb lattice 
is obviously larger than the one on the square lattice 
with J cl = 0.56A < J < J c2 = 0.62AJ£ We suspect 
this difference may result from unique linear DOS on the 
honeycomb lattice at half-filling, which lowers the energy 
penalty of coexistence. 

IV. EXTENSIONS AND DISCUSSIONS 

Moreover, for physical observable, it is expect that 
in the proposed coexistence regime the dynamical mag- 
netic structure factor, which can be detected in inelastic 
neutron scattering experiments, should have two contri- 
butions. One is from a momentum-independent single 
site slow component due to the localized Kondo-type ex- 
citations, and the other is contributed from a strongly 
momentum-dependent inter-site fast component, reflect- 
ing the magnetic interactions^ 

We noted that, it is also interesting to study a more 
simplified model only including degrees of freedom for 
spins, namely, the Kondo necklace model on the honey- 
comb latticed However, the present mean-field decou- 
pling approach is not cheap in such a spin model but the 
so-called bond-operator representation could be a use- 
ful tool in this case4£ Moreover, if a spin-orbit coupling 
is introduced in our present model as the case in Ref. 
[36| , one could study the interplay among the antifer- 



romagnetic order, Kondo screening and the symmetry- 
protected topological states (such as the 2D topological 
insulator). But, it seems that the topological insulating 
state found in Ref. [36[ might be unstable to the con- 
ventional antiferromagnetic order in terms of the present 
mean-field decoupling treatment^ Therefore, more so- 
phisticated treatments on this issue is eagerly desirable. 



V. CONCLUSION 

In summary, we have obtained the ground-state phase 
diagram of the Kondo lattice model on the honeycomb 
lattice at half-filling using extended mean-field decou- 
pling. By treating magnetic interaction and Kondo 
screening on an equal footing, it is found that besides 
a trivial discontinuous first-order quantum phase transi- 
tion between well-defined Kondo insulator and antiferro- 
magnetic insulating state, there can exist a wide coexis- 
tence regime with both Kondo screening and antiferro- 
magnetic orders in the intermediate coupling regime. It 
is expected such a coexistence regime could be realized 
by experiments of ultra-cold atoms on the honeycomb 
optical lattices and may be found by quantum Monte 
Carlo simulation in future since it is free of the 'minus 
sign'i 16 ' 17 i 42 i 43 We hope the present work may be helpful 
for further studies on the interplay between conduction 
electrons and the densely localized spins for the honey- 
comb lattice. 
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